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Abstract — Subsystem codes protect quantum information by 
encoding it in a tensor factor of a subspace of the physical 
state space. Subsystem codes generalize all major quantum error 
protection schemes, and therefore are especially versatile. This 
paper introduces numerous constructions of subsystem codes. It 
is shown how one can derive subsystem codes from classical 
cyclic codes. Methods to trade the dimensions of subsystem 
and co-subystem are introduced that maintain or improve the 
minimum distance. As a consequence, many optimal subsystem 
codes are obtained. Furthermore, it is shown how given subsystem 
codes can be extended, shortened, or combined to yield new 
subsystem codes. These subsystem code constructions are used 
to derive tables of upper and lower bounds on the subsystem 
code parameters. 



I. Introduction 

Quantum information processing as a growing exciting field 
has attracted researchers from different disciplines. It utilizes 
the laws of quantum mechanical operations to perform expo- 
nentially speedy computations. In an open system, one might 
wonder how to perform such computations in the presence 
of decoherence and noise that disturb quantum states storing 
quantum information. Ultimately, the goals of quantum error- 
correcting codes are to protect quantum states and to allow 
recovery of quantum information processed in computational 
operations of a quantum computer Henceforth, one seeks to 
design good quantum codes that can be efficiently utilized for 
these goals. 

A well-known approach to derive quantum error-correcting 
codes from self-orthogonal (or dual-containing) classical codes 
is called stabilizer codes, which were introduced a decade ago. 
The stabilizer codes inherit some properties of Clifford group 
theory, i.e., they are stabilized by abelian finite groups. In the 
seminal paper by Calderbank at. et [7], [20], [22], various 
methods of stabilizer code constructions are given, along with 
their propagation rules and tables of upper bounds on their 
parameters. In a similar tactic, we also present subsystem 
code structures by establishing several methods to derive them 
easily from classical codes. Subsystem codes inherit their 
name from the fact that the quantum codes are decomposed 
into two systems as explained in Section II. The classes of 
subsystem codes that we will derive are superior because 
they can be encoded and decoded using linear shirt-register 
operations. In addition, some of these classes turned out to be 
optimal and MDS codes. 



Subsystem codes as we prefer to call them were mentioned 
in the unpublished work by Knill [14], [15], in which he 
attempted to generalize the theory of quantum error-correcting 
codes into subsystem codes. Such codes with their stabiUzer 
formalism were reintroduced recently [6], [12], [16], [17], 
[ ]. An ((n, iC, R, d))q subsystem code is a .fsT i?-dimensional 
subspace Q of C that is decomposed into a tensor product 
Q = A® B of a. if-dimensional vector space A and an R- 
dimensional vector space B such that all errors of weight less 
than d can be detected by A. The vector spaces A and B are 
respectively called the subsystem A and the co-subsystem B. 
For some background on subsystem codes see the next section. 

This paper is structured as follows. In section U, we present 
a brief background on subsystem code structures and present 
the Euclidean and Hermitian constructions. In section UI, 
we derive cyclic subsystem codes and provide two generic 
methods of their constructions from classical cyclic codes. 
Consequently in section IV, we construct families of subsys- 
tem BCH and RS codes from classical BCH and RS over 
and defined using their defining sets. In Sections V,VI,Vn, 
we establish various methods of subsystem code constructions 
by extending and shortening the code lengths and combining 
pairs of known codes, in addition, tables of upper bounds on 
subsystem code parameters are given. Finally, the paper is 
concluded with a discussion and future research directions in 
section VIII. 

Notation. If is a set, then \S\ denotes the cardinality of 
the set S. Let g be a power of a prime integer p. We denote 
by Fg the finite field with q elements. We use the notation 
{x\y) = (xi, . . . , Xn\yl^ . . . , Un) to dcnotc the concatenation 
of two vectors x and y in F^. The symplectic weight of 
{x\y) e F^" is defined as 

swt{x\y) = {{x,,y^) ^ (0,0) \l<i< n}. 

We define swt(X) = min{swt(a;) \ x ^ X,x ^ 0} for any 
nonempty subset X ^ {0} of F^". 

The trace-symplectic product of two vectors u — (a|6) and 
V = {a'\b') in F^" is defined as 

{u\v)s = iyq/p{a - b- a - 6'), 

where x ■ y denotes the dot product and tr^/p denotes the 
trace from F^ to the subfield Fp. The trace-symplectic dual 
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of a code C C F^" is defined as 

C^^ = {w e F^" I {v\w)s = for all w G C}. 

We define the Euclidean inner product {x\y) = X^ILi ^iVi ^^'^ 
the Euclidean dual of C C F^"^ as 

^{xe F;' I {x\y) = for all y e C}. 

We also define the Hermitian inner product for vectors x, y 
in F^2 as {x\y)h ~ X)r=i y« ^'^'^ '■^^ Hermitian dual of 
C C F;\ as 

C^'^ ={xe F^. I {x\y)h = for all y 6 C}. 

II. Background on Subsystem Codes 

In this section we give a quick overview of subsystem codes. 
We assume that the reader is familiar the theory of stabilizer 
codes over finite fields, see [7], [II], [20] and the references 
therein. 



A. Errors 

Let Fg denote a finite field with q elements of characteris- 
tic p. Let {|a-) I x G Fg} be a fixed orthonormal basis of C 
with respect to the standard hermitian inner product, called 
the computational basis. For a, 6 G Fg, we define the unitary 
operators X{a) and Z{h) on by 

X{a)\x) = |.T + a), Z{b)\x) = 

where uj = exp(27ri/p) is a primitive pth root of unity 
and tr is the trace operation from F^ to Fp. The set E = 
{X{a)Z{b)\a,h G F^} forms an orthogonal basis of the 
operators acting on with respect to the trace inner product, 
called the error basis. 

The state space of n quantum digits (or qudits) is given 
by C«" = C? ® C? ® • • • (g) C. An error basis E on C«" 
is obtained by tensoring n operators in E; more explicitly, 
E = {X(a)Z(b) I a,b G F^}, where 

X{sl) = X{ai)®---®X{an), 
Z{h) = Z{bi) ® • • • ® Z{bn) 

for a = (ai, . . . , a„) G F^ and b = (hi, . . . , 6„) G F^. The 
set E is not closed under multiplication, whence it is not a 
group. The group G generated by E is given by 

G = {o.'^E = uj'X{a)Z{h) I a, b G F^ , c G Fp}, 

and G is called the error group of C . The error group is 
an extraspecial p-group. The weight of an error in G is given 
by the number of nonidentity tensor components; hence, the 
weight of uj'^X{a)Z{h) is given by the symplectic weight 
swt(a|b). 



B. Subsystem Codes 

An ((n, R, d))q subsystem code is a subspace Q = A(E)B 
of C that is decomposed into a tensor product of two vector 
spaces A and B of dimension dim A = K and dimi? = R 
such that all errors in G of weight less than d can be detected 
by A. We call A the subsystem and B the co-subsystem. The 
information is exclusively encoded in the subsystem A. This 
yields the attractive feature that errors affecting co-subsystem 
B alone can be ignored. 

A particularly fruitful way to construct subsystem codes 
proceeds by choosing a normal subgroup N of the error 
group G, and this choice determines the dimensions of sub- 
system and co-subsystem as well as the error detection and 
correction capabilities of the subsystem code, see [12]. One 
can relate the normal subgroup to a classical code, namely 
modulo the intersection of A^ with the center Z{G) of 
G yields the classical code X = N/{N n Z{G)). This 
generalizes the familiar case of stabilizer codes, where A^ is 
an abelian normal subgroup. It is remarkable that in the case 
of subsystem codes any classical additive code X can occur 
It is most convenient that one can also start with any classical 
additive code and obtain a subsystem code, as is detailed in 
the following theorem from [12]: 

Theorem 1. Let C be a classical additive subcode of F^" 
such that C ^ {0} and let D denote its subcode D = CnC^°. 
If X = |C| and y = \D\, then there exists a subsystem code 
Q ^ A® B such that 

i) dimA = q"/{xy)^/'^, 

ii) diuiB = {x/yy/^. 

The minimum distance of subsystem A is given by 

(a) d ^ swt((C + C-L^ -C)= swt{D^^^ - C) if D^^^ ^ G; 

(b) d ^ swtp^-) if D^^ = C. 

Thus, the subsystem A can detect all errors in E of weight less 
than d, and can correct all errors in E of weight < [(d— 1)/2J. 

Proof: See [12, Theorem 5]. ■ 
A subsystem code that is derived with the help of the 
previous theorem is called a Clifford subsystem code. We 
will assume throughout this paper that all subsystem codes 
are Clifford subsystem codes. In particular, this means that 
the existence of an {{n, K, R,d))g subsystem code implies 
the existence of an additive code C < F^" with subcode 
D = C n C^" such that |C| = q"R/K, \D\ = q"/{KR), 
and d = swt{D^'' - C). 

A subsystem code derived from an additive classical code C 
is called pure to d' if there is no element of symplectic weight 
less than d' in C. A subsystem code is called pure if it is pure 
to the minimum distance d. We require that an ((n, 1, R, d))q 
subsystem code must be pure. 

We also use the bracket notation [[n, fc,r, d]]g to write the 
parameters of an {{n,q^ ,q^ ,d))q subsystem code in simpler 
form. Some authors say that an [[n, fc, r, djjg subsystem code 
has r gauge qudits, but this terminology is slightly confusing, 
as the co-subsystem typically does not correspond to a state 
space of r qudits except perhaps in trivial cases. We will avoid 
this misleading terminology. An {{n,K,l,d))q subsystem 
code is also an {{n, K, d))q stabilizer code and vice versa. 
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Subsystem codes can be constructed from the classical 
codes over Fg and F^2 . We recall the Euclidean and Hermitian 
constructions from [j], which are easy consequences of the 
previous theorem. 

Lemma 2 (Euclidean Construction). If C is a k' -dimensional 
Fq-linear code of length n that has a k" -dimensional subcode 
D = C n C'^ and k' + k" < n, then there exists an 

[[n, n-{k' + k"), k' - k", wt{D^ \ C)]]g 

subsystem code. 

Lemma 3 (Hermitian Construction). If C is a k' -dimensional 
Fg2-linear code of length n that has a k" -dimensional subcode 
D = C D C^'' and k' + k" < n, then there exists an 

[[n, n-{k' + k"),k' - k", wt{D^'^ \ C)]]g 

subsystem code. 

III. Cyclic Subsystem Codes 

In this section we shall derive subsystem codes from 
classical cyclic codes. We first recall some definitions before 
embarking on the construction of subsystem codes. For further 
details concerning cyclic codes see for instance [10] and [18]. 

Let n be a positive integer and a finite field with q 
elements such that gcd(n, q) = 1. Recall that a linear code 
C C Fg is called cyclic if and only if (cg, . . . , c„_i) in C 
implies that (c„_i, cg, . . . , c„_2) in C. 

For g{x) in Fg[x], we write {g{x)) to denote the principal 
ideal generated by g{x) in [x]. Let tt denote the vector space 
isomorphism tt: F^ ^ i?„ = Fg[x]/{x" — 1) given by 

7r((co, . . . , c„_i)) = Co + cia; H h c„_ia;"~^ + (x" - 1). 

A cyclic code C C F^'^ is mapped to a principal ideal tt{C) of 
the ring i?„. For a cyclic code C, the unique monic polynomial 
g{x) in Fg[x] of the least degree such that {g{x)) = 7r(C) is 
called the generator polynomial of C. If C C F^ is a cyclic 
code with generator polynomial g{x), then 

dimFg C = n — degg{x). 

Since gcd(n, q) = 1, there exists a primitive n* root of 
unity a over F,; that is, Fq\a\ is the splitting field of the 
polynomial x" — 1 over Fg. Let us henceforth fix this primitive 
n* primitive root of unity a. Since the generator polynomial 
g{x) of a cyclic code C C F^ is of minimal degree, it follows 
that g{x) divides the polynomial x" — 1 in Fg[a;]. Therefore, 
the generator polynomial g{x) of a cyclic code C C F^ can 
be uniquely specified in terms of a subset T of {0, . . . , n — 1} 
such that 

g{x) = \{{x-a'). 

The set T is called the defining set of the cyclic code C (with 
respect to the primitive n* root of unity a). Since g{x) is a 
polynomial in Fg[a;], a defining set is the union of cyclotomic 
cosets Cx, where 

Cx = {xq"" mod n I i e Z, i > 0}, < a; < n 



The following lemma recalls some well-known and easily 
proved facts about defining sets (see e.g. [10]). 

Lemma 4. Let Ci be a cyclic code of length n over Fg with 
defining set a Ti for i = 1,2. Let N = {0, 1, . . . , n — 1} and 
Ti ~ {at mod n\t E T} for some integer a. Then 

i) Ci n C'2 has defining set Ti U T2. 

ii) Ci + C2 has defining set Ti n T2. 
Hi) C*! C C2 if and only if T2 C Ti. 

iv) has defining set N \ T^^. 

v) Cj^'' has defining set N \ T^"^ provided that q ~ for 
some positive integer r. 

Notation. Throughout this section, we denote by N the set 
A^ = {0,...,n— 1}. The cyclotomic coset of x will be denoted 
by Ca;. If T is a defining set of a cyclic code of length n and 
a is an integer, then we denote henceforth by the set 

= {at mod n\t(E T}, 

as in the previous lemma. We use a superscript, since this no- 
tation will be frequently used in set differences, and arguably 
TV \ is more readable than N\-qT. 

Now, we shall give a general construction for subsystem 
cyclic codes. We say that a code C is self-orthogonal if and 
only if C C C^. We show that if a classical cyclic code is 
self-orthogonal, then one can easily construct cyclic subsystem 
codes. 

Proposition 5. Let D be a k-dimensional self- orthogonal 
cyclic code of length n over Fg. Let T]j and Tjj± respectively 
denote the defining sets of D and D^. If T is a subset of 
To \ To± that is the union of cyclotomic cosets, then one can 
define a cyclic code C of length n over Fg by the defining 
set Tc ^TdMTU T-1). If r = \T U T"!] is in the range 
< r < n — 2k, and d = minwt(D^ \ C), then there exists 
a subsystem code with parameters [[n, n — 2fc — r, r, d]]g. 

Proof: Since 13 is a self-orthogonal cyclic code, we have 
D C D^, whence T]j± C by Lemma 4 iii). Observe that 
if s is an element of the set S = Td\Td± =Td\ {N\T^^), 
then — s is an element of S as well. In particular, is a 
subset of Td\Tjj±. 

By definition, the cyclic code C has the defining set Tc = 
Td \ (T U T^i); thus, the dual code has the defining set 

Tc± TJ^ Ta± U (T U T-i). 

Furthermore, we have 

Tc U Tc± = [To \ (T U T-i)) U (T^^ U T U T'^) = T^; 

therefore, C n = by Lemma 4 i). 

Since n-k ^ \Td \ and r = |ruT"i|, we have dimp, D = 
n - |Td| = k and dimp, C ^ n - \Tc\ = k + r. Thus, 
by Lemma 2 there exists an F, -linear subsystem code with 
parameters [[n, k, p, where 

i) K = dim D-^ — dim C = rt — fc — (fc + r) = n — 2k — r, 

ii) p — dim C — dim D = k + r — k = r, 

iii) d = minwt(D-L \ C), 

as claimed. ■ 
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We can also derive subsystem codes from cyclic codes over 
Fq2 by using cyclic codes that are self-orthogonal with respect 
to the Hermitian inner product. 

Proposition 6. Let D be a cyclic code of length n over Fg2 
such that D C Let Tjj and Tj~i±,^ respectively be the 

defining set of D and D^*^. If T is a subset of To \ Tjj±^ 
that is the union of cyclotomic cosets, then one can define a 
cyclic code C of length n over Fg2 with defining set Tq = 
Td\{TU T-1). If n - k = \Td\ and r = \T U T"?! with 
< r < n — 2k, and d = wt(£)-^'' \ C), then there exists an 
[[n,n — 2k — r,r,d]]q subsystem code. 

Proof: Since D C D^'\ their defining sets satisfy 
Tpi^ C To by Lemma 4 iii). If s is an element of To\To±h, 
then one easily verifies that —qs (mod n) is an element of 

To \ Tjj±,^ . 

Let N = {0, 1, . . . , n — 1}. Since the cyclic code C has the 
defining set Tc ^ To \ [T U T"'), its dual code C^" has 
the defining set Tc±,, = N XT^'' = To±,, U (T U T""?). We 
notice that 

Tc U Tc±, = {To \ (T U T-")) U (T^^, U T U T"') = To; 

thus, C n C^** = Dhy Lemma 4 i). 

Since n-k= \To\ and r^\T\J T"«|, we have dimD = 
n — \To\ — k and dimC = n — \Tc\ ~ k + r. Thus, by 
Lemma 3 there exists an [[n, k, p, d]]q subsystem code with 

i) K = dim £)-'-'• — dim C = {n — k) — {k + r) ~ ri — 2fc — r, 

ii) p = dim C — dim D^k + r — k — r, 
iii) d = minwt(L»-L'' \ C), 

as claimed. ■ 
We include an example to illustrate the construction given 
in the previous proposition. 

Example 7. Consider the narrow-sense BCH code D^i- of 
length n = 31 over F4 with designed distance 5. The defining 
set To±H of D-^'- is given by To±h = Ci U C2 U C3 U C4 = 
Ci U C3, where the cyclotomic cosets of 1 and 3 are given by 

Ci = {1,4, 16,2,8} and C3 = {3, 12, 17, 6, 24}. 

If N = {0, 1, . . . , 30}, then the defining set of the dual code D 
is given byTo=N\ {C15 U C7) = Co U Ci U C3 U C5 U Cn . 
Therefore, D C D^'\ dimD-L" = 21 and dimD = 10. // 
we choose T = C5, then T^^ = Cn, whence the defining 
set Tc of the code C is given by Tc = To \ (C5 U Cn) = 
Co UCi UC3. It follows that dimC 20 and dimC-^" = 11. 
Therefore, the construction of the previous proposition yields 
a BCH subsystem code with parameters [[31, 1, 10, > 5]]2. 

The general principle behind the previous example yields 
the following simple recipe for the construction of subsystem 
codes: Choose a cyclic code (such as a BCH or Reed-Solomon 
code) with known lower bound 5 on the minimum distance 
that contains its (hermitian) dual code, and use Proposition 5 
(or Proposition 6) to derive subsystem codes. This approach 
allows one to control the minimum distance d of the subsystem 
code, since d > 5 is guaranteed. Another advantage is that 
one can exploit the cyclic structure in encoding and decoding 
algorithms. 



TABLE I 

SUBSYSTEM BCH CODES THAT ARE DERIVED USING THE EUCLIDEAN 
CONSTRUCTION 



OUUSyalClll ^^UU-C 




L'CftlgllCU. 






Hi (itfinff 


[[10, 4t, v5, v3Jj2 


[iO, / , OJ2 


4 


[[±0, U, -L, OJJ2 


1 c; c; 7 L, 

[± J, 0, ( J2 




[[31,10,1,5]]2 


[31,11,11]2 


8 


20 1 So 


[O J- , U, -LiJJ 2 


12 


[[DO, D, Zl, /JJ2 


[fiS "^Q Qlo 
[UO , oy , yj 2 




[[uo, u, itj, 1 JJ 2 


[UO , OU, J- ±j 2 


10 


[[UO, 0, ' JJ i 


fiS SO 

IVJiJ, OVJ, J.tJJ^ 


12 


\\(\'\ IS '\ 7lln 
[[UO, ±0, 0, /JJ2 


ffi'^ 94 1 '^lo 

[UO , Zt:, J-OJ 2 


14 


[[UO, ou, 0, iJJJ 2 


[fis IS 9110 

[UO , J-O, ^ ±J 2 


16 


[[UO, Oii, ±, "JJJ 2 


ffi'^ Ifi 9'^lo 

[UO , ±U, ZOJ 2 


22 


[[UO, ^^4:, ±, OJJ 2 


fiS in 97 


24 


[[63, 50, 1, 3]]2 


[63, 7, 31] 2 


28 


[[15,2,5,3]]4 


[15,9,5]4 


4 


[[15,2,3,3]]4 


[15,8,6]4 


6 


[[15,4,1,3]]4 


[15,6,7]4 


7 


[[15,8,1,3]]4 


[l5,4, 10]4 


8 


[[31,10, 1,5]]4 


[31,11,11]4 


8 


[[31,20, 1,3]]4 


[31,6, 15]4 


12 


[[63,12,9,7]]4 


[63,30, 15]4 


15 


[[63,18,9,7]]4 


[63,27, 21]4 


16 


[[63,18,7, 7]]4 


[63,26, 22]4 


22 



* punctured code 
+ Extended code 



For example, if we start with primitive, narrow-sense BCH 
codes, then Proposition 5 yields the following family of 
subsystem codes: 

Corollary 8. Consider a primitive, narrow-sense BCH code 
of length n — q™ — 1 with m > 2 over Fq with designed 
distance 6 in the range 

2<S < gr™/2l _ 1 _ _ 2)[m is odd]. (1) 

If T is a subset of N \ { [J^a^\{Ca U C_a)) that is a union of 
cyclotomic cosets and r = \T U T~^\ with < r < n — 2k, 
where k ~ m\{S — 1)(1 — l/q)~\, then there exists an 

[[g™ ~ 1, - 1 - 2m\i6 - 1)(1 - 1/q)] - r, r, > 5]]q 

subsystem code. 

Proof: By [ Theorem 2], a primitive, narrow-sense BCH 
code D-^ with designed distance 6 in the range (1) satisfies 
D C D-^. By [5, Theorem 7], the dimension of is given 
by dimZ?^ = g'" - 1 - m.\{S - 1)(1 - l/q)] = n - k, 
whence k ~ dim D. Let To and To± respectively denote the 
defining sets of D and D^. It follows from the definitions that 
To± = Ua=i and that T is a subset of 

N \ {To^ UT^l) = {N\ T^i) \ To^ = To \ To^. 

If Tc ^Td\{TUT-^) denotes the defining set of a cycUc 
code C, then dim C = k + r. By Proposition 5, there exists 
an [[n, n — 2fc — r, r, > 6]]q subsystem code, which proves the 
claim. ■ 
Similarly, we can obtain a hermitian variation of the pre- 
ceding corollary with the help of Proposition 6. 
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TABLE II 

Subsystem BCH codes that are derived with the help of the 
hermitian construction 



Subsystem Code 


Parent BCH 
Code C 


Designed 
distance 


1[14,1,3,4JJ2 
[[15,1,2,5]]2 
[[15,5,2,3]]2 
[[16,5,2,3]]2 


[14,8,5]22 
[15,8,6]22 
[15,6,7]22 
[16,6,7]22 


6* 
6 
7 

7+ 


1[17,8,1,4JJ2 


[17,5,9)22 


4 


1[21,6,3,3JJ2 
[[21,7,2,3]]2 


[21,9,7]]22 
[21,8,9]22 


6 
8 


[[31,10,1,5||2 
[[31,20,1,3]]2 
[[32,10,1,5]]2 
[[32,20,1,3]]2 


[31,11,11J22 

[31,6,15]22 
[32,11,11)22 
[32,6, 15]22 


8 

12 

8+ 
12+ 


[[25,12,3,3] 
[[26,6,2,5]] 
[[26,12,2,4] 
[[26,13,1,4] 


3 
3 
3 
3 


[25, 8,12)32 
[26, 11,8)32 

[26.8. 13) 32 

[26.7. 14) 32 


9* 

8 

9 
14 


[[80,1,17, 20 
[[80,5, 17, 17 


]3 
]3 


[80,48,21]32 
[80,46,22]32 


21 
22 



* punctured code 
+ Extended code 



Corollary 9. Consider a primitive, narrow-sense BCH code 
of length n ~ g^™ — 1 with to ^ 2 over Fg with designed 
distance 5 in the range 

2 < (5 < - 1 (2) 

If T is a subset of the set N \ ^lJa=i(^a C'-qa)^ that is a 
union of cyclotomic cosets and r = |T U T^'^l with < r < 
n — 2k, where k = m\(S — 1)(1 — then there exists a 

subsystem code. 

Proof: The proof is similar to the proof of the previous 
corollary, and is a consequence of [5, Theorems 4 and 7] and 
Proposition 6. ■ 

It is straightforward to generalize the previous two corollar- 
ies to the case of non-primitive BCH codes using the results 
given in [4], [2]. 

One of the disadvantages of the cyclic constructions is that 
the parameter r is restricted to values dictated by the possible 
cardinaUties of the sets T U T"! (or T U T'l), where T 
is confined to be a union of cyclotomic cosets. In the next 
section, we will see how one can overcome this limitation. 

We conclude this section by giving some examples of the 
parameters of subsystem BCH codes in Tables I and II. 

IV. Trading Dimensions of subsystem and 

CO-SUBSYSTEM CODES 

In this section we show how one can trade the dimensions 
of subsystem and co-subsystem to obtain new codes from a 
given subsystem or stabilizer code. The results are obtained by 
exploiting the symplectic geometry of the space. A remarkable 
consequence is that nearly any stabilizer code yields a series 
of subsystem codes. 

Our first result shows that one can decrease the dimension 
of the subsystem and increase at the same time the dimension 



of the co-subsystem while keeping or increasing the minimum 
distance of the subsystem code. 

Theorem 10. Let q be a power of a prime p. If there exists an 
{{n, K, R, d))q subsystem code with K > p that is pure to d' , 
then there exists an [in, K/p,pR, > d))q subsystem code that 
is pure to min{d, d'}. If a pure {{n^p, R, d))q subsystem code 
exists, then there exists a {{n,l,pR,d))q subsystem code. 

Proof: By definition, an {{n, K, R, d))q Clifford subsys- 
tem code is associated with a classical additive code C C F^" 
and its subcode D = C n such that x = |C|, y = \D\, 
K = g"/(.Ty)i/2, R = {x/y^/^, and d = swt{D^^ - C) if 
C ^ D^=, otherwise d = swt{D^^) if = C. 

We have q = for some positive integer to. Since K and 
R are positive integers, we have x = p^+2'' and y ~ p'^ for 
some integers r > 1, and s > 0. There exists an Fp-basis of 
C of the form 

C = Spailp^jzi, . . . , Zg, Xg+l, Zs+li ■ • ■ 1 ^s+r, Zs+r} 

that can be extended to a symplectic basis 

{xi,Zi,...,Xnm,Zn,n} of F^", that is, (xfc | Xe)s = 0, 

{zk I Zi)s = 0, (xfc I Zi)s = 5k,i for all 1 < /c,^ < nm, 
see [8, Theorem 8. 10. 1]. 
Define an additive code 

Cm ~ Spailp^jzi, . . . ,Zs, Xs+1: ■ ■ • i Xs+r+l: ^s+r+l}- 

It follows that 
and 

D = Cm n C^" = spanp^{zi, . . . , z^}. 

By definition, the code C is a subset of Cm- 

The subsystem code defined by Cm has the parameters 

in,Km,Rm,dm), whcrC Km - q"/{p-+^r+2pS^l/2 ^ ^/^ 

and Rm ~ = pR. For the claims concerning 

minimum distance and purity, we distinguish two cases: 

(a) If Cm then K > p and d,„ swt(i:>^' - C,„) > 
swt{D^= -C) = d. Since by hypothesis swt(£>^' -C) = 
d and swt(C) > d' , and £> C C C C„ C D^^ by 
construction, we have swt(Cm) > minjii, d'}; thus, the 
subsystem code is pure to min{d, d'}. 

(b) If Cm = -D^% then Km = 1 = K/p, that is, K = p; it 
follows from the assumed purity that d = swt(Z?^= — C) = 

SWt(Z?-'-=) = dm- 

This proves the claim. ■ 
For Fq-linear subsystem codes there exists a variation of 
the previous theorem which asserts that one can construct the 
resulting subsystem code such that it is again Fg-linear. 

Theorem 11. Let q be a power of a prime p. If there exists an 
Fq-linear [[n, k, r, d]]q subsystem code with k > 1 that is pure 
to d', then there exists an Fq-linear [[n,k — l,r + 1,> d\]q 
subsystem code that is pure to minjd, d'}. If a pure Fq-linear 
[[n,l,r,d]]q subsystem code exists, then there exists an Fq- 
linear [[n, 0, r + 1, d]]q subsystem code. 

Proof: The proof is analogous to the proof of the previous 
theorem, except that F^-bases are used instead of Fp-bases. 
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There exists a partial converse of Theorem 10, namely if 
the subsystem code is pure, then it is possible to increase the 
dimension of the subsystem and decrease the dimension of the 
co-subsystem while maintaining the same minimum distance. 

Theorem 12. Let q be a power of a prime p. If there exists 
a pure {{n, K, R, d))q subsystem code with R > 1, then there 
exists a pure ({n,pK, R/p,d))q subsystem code. 

Proof: Suppose that the {{n, K, R, d))q Clifford subsys- 
tem code is associated with a classical additive code 

Cm = SpailFpl^l, . . . , Zs, Xs+1, Zs+1, . . . , Xs+r+1, ^s+r+l}- 

Let D = CmD C,i^. We have x = |C„| = y = 

\D\ = p*, hence K = (7"/p'"+'* and R = Furthermore, 

The code 

C = Spanp^{2;i, . . . , Z^, Xs+l, Zs+l, . . . ,Xs+r, Z^+r} 

has the subcode D ^ C n C^=. Since |C| = |Cm|/p^, the 
parameters of the Clifford subsystem code associated with C 
are {{n,pK, R/p, d'))q. Since C C Cm, the minimum distance 
d' satisfies 

d' = swt(L>-L-' -C) < swt(7:>-L'- - C,„) = swt(i:>^-') = d. 

On the other hand, d' = swt(-D-^= - C) > swt{D-^-) = d, 
whence d = d'. Furthermore, the resulting code is pure since 
d = swt{D^^ ) = swt(L>-L= - C). m 
Replacing Fp-bases by F^-bases in the proof of the previous 
theorem yields the following variation of the previous theorem 
for Fg-linear subsystem codes. 

Theorem 13. Let q be a power of a prime p. If there exists a 
pure Fq-linear [[n^ k,r, d]]q subsystem code with r > 0, then 
there exists a pure Fq-linear [[n, A; + 1, r — 1, d]]q subsystem 
code. 

The purity hypothesis in Theorems 12 and 13 is essential, 
as the next remark shows. 

Remark 14. The Bacon-Shor code is an impure [[9, 1, 4, 3]]2 
subsystem code. However, there does not exist any [[9, 5,3]]2 
stabilizer code. Thus, in general one cannot omit the purity 
assumption from Theorems 12 and 13. 

An [[n, fc, d\]q stabilizer code can also be regarded as an 
[[n, fc, 0, d\]q subsystem code. We record this important special 
case of the previous theorems in the next corollary. 

Corollary 15. If there exists an (Fq-linear) [[n, fc,d]]q stabi- 
lizer code that is pure to d', then there exists for all r in 
the range < r < k an (Fq-linear) [[n,k — r,r,> d]]q 
subsystem code that is pure to Tam{d,d'} . If a pure (Fq- 
linear) [[n, k,r, d]]q subsystem code exists, then a pure (Fq- 
linear) [[n, fc + r, d]]q stabilizer code exists. 

We have shown in [4], [5] that a (primitive or non-primitive) 
narrow sense BCH code of length n over Fq contains its dual 
code if the designed distance 6 is in the range 

2<s< <5„ax = -^{q^"'^'^ -i-{q- 2)[m odd]). 



For simplicity, we will proceed our work for primitive nar- 
row sense BCH codes, however, the generalization for non- 
primitive BCH codes is a straightforward. 

Corollary 16. If q is power of a prime, m is a positive integer, 
and 2 < S < (/T™/^! — 1 — (q — 2)[m odd ]. Then there exists a 
subsystem BCH code with parameters [[g™ — 1, n — 2m\{5 — 

1) (l-l/g)l -r,r, > 5]], where < r < n-2m\{5 

Proof We know that if 2 < J < qTm/sl - l - - 

2) [m odd ], then there exists a stabilizer code with parameters 
[[q™ - 1, 71 - 2m\{5 - 1)(1 - l/q)l , > 5]]q. Let r be an integer 
in the range < r < n — 2m\{5 — 1)(1 — I/'?)]- From [1, 
Theorem 2], then there must exist a subsystem BCH code with 
parameters [[g™ — 1, n — 2m\{5 — 1)(1 — 1/<7)1 — > <5]]q- 

■ 

We can also construct subsystem BCH codes from stabilizer 
codes using the Hermitian constructions. 

Corollary 17. If q is a power of a prime, m is a positive 
integer, and S is an integer in the range 2 < S < (5max = 
qm+lm even]_-[^_^g2 2)[77i evcn], then there exists a subsystem 
code Q with parameters 

[fe"" - 1, - 1 - 2m\{5 - 1)(1 - l/q^)} - r, r, > 5% 
that is pure up to 5, where < ?' < q^"^ — 1 — 2m [(5 — 1)(1 — 

Proof If 2 < (5 < (5„iax = ""'"l - I - {q^ - 

2)[m even], then exists a classical BCH code with parameters 
[g'"-l,g'"-l-m[(J-l)(l-l/g)],> S\q which contains 
its dual code. From [I, Theorem 2], then there must exist a 
subsystem code with the given parameters. ■ 

V. MDS Subsystem Codes 

Recall that an [[ri, fc, r, d]]g subsystem code derived from 
an Fq-linear classical code C < Fq"' satisfies the Singleton 
bound k + r < n~~2d+2, see [13, Theorem 3.6]. A subsystem 
code attaining the Singleton bound with equality is called an 
MDS subsystem code. 

An important consequence of the previous theorems is the 
following simple observation which yields an easy construc- 
tion of subsystem codes that are optimal among the F^-linear 
Clifford subsystem codes. 

Theorem 18. If there exists an Fq-linear [[n, k, d]]q MDS sta- 
bilizer code, then there exists a pure Fq-linear [[n, k — r, r, d]]q 
MDS subsystem code for all r in the range < r < k. 

Proof: An MDS stabilizer code must be pure, see [20, 
Theorem 2] or [11, Corollary 60]. By Corollary 15, a pure 
Fg-linear [[n, k, d]]q stabiUzer code implies the existence of an 
Fg-linear [[n,k — r,r,dr > d]]q subsystem code that is pure 
to d for any r in the range < r < fc. Since the stabilizer code 
is MDS, we have k = n — 2d+2. By the Singleton bound, the 
parameters of the resulting Fg-linear [[n,n — 2d+2 — r, r, dr]]q 
subsystem codes must satisfy {n — 2d-\-2 — r)+r< n~2d,. + 2, 
which shows that the minimum distance d^ = d, as claimed. 
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Remark 19. We conjecture that Fq-linear MDS subsystem 
codes are actually optimal among all subsystem codes, but a 
proof that the Singleton bound holds for general subsystem 
codes remains elusive. 

In the next corollary, we give a few examples of MDS sub- 
system codes that can be obtained from Theorem 18. These are 
the first families of MDS subsystem codes (though sporadic 
examples of MDS subsystem codes have been established 
before, see e.g. [3], [6]). 

Corollary 20. i) An Fq-linear pure [[n, n—2d+2~r, r, d]]q 
MDS subsystem code exists for all n, d, and r such that 
3 < n < q, 1 < d < n/2 + 1, and < r < n - 2d + 1. 

ii) An Fq-linear pure [[{h'+l)q, {v+l)q—2i'—2—r,r,v+2]]q 
MDS subsystem code exists for all v and r such that < 
V <q-2 and'^i <r <(y + \)q-2v - i. 

Hi) An Fq-linear pure [[q — 1, q — 1 — 26 — r,r, 6 + l]]q MDS 
subsystem code exists for all 6 and r such that < d < 
{q-l)/2andO<r<q-2S-l. 

iv) An Fq-linear pure [[q, q ~ 25 - 2 ~ r' ,r' ,5 + 2\]q MDS 
subsystem code exists for all Q < 5 < (q — l)/2 and 
< r' < g - 2(5 - 2. 



TABLE III 
Optimal pure subsystem codes 



v) An Fq-linear pure [[q^ — — 2(5 — 1 — 



1]]. 



MDS subsystem code exists for all 6 and r in the range 
< (5 < g - 1 fl«<i < r < g2 _ 2(5 - 1. 
vi) An Fq-linear pure [[q'^, q'^ -26-2- r' , r' , S + 2]], MDS 
subsystem code exists for all 5 and r' in the range < 
5 < q-1 and Q <r' < q^ -25 -2. 

Proof: i) By [9, Theorem 14], there exist Fq-linear 
[[n, 71 — 2(i + 2, d\]q stabilizer codes for all n and d such that 
i < n < q and 1 < d < n/2 + 1. The claim follows from 
Theorem 18. 

ii) By [21, Theorem 5], there exist a [[(z^ + l)q, [v + \)q — 
21/ — 2,1/ + 2]]q stabilizer code. In this case, the code is 
derived from an Fg2 -linear code X of length n over Fq2 
such that X C X-^'\ The claim follows from Lemma 29 and 
Theorem 18. 

iii) ,iv) There exist Fq-linear stabilizer codes with parameters 

[[q-l,q-2S -l,S+l]]q and [[q,q - 2S - 2,6 + 2]]q for 
< 5 < (q — 1) /2, see [ , Theorem 9]. Theorem 18 yields 
the claim. 

v),vi) There exist Fg-linear stabilizer codes with parameters 
[[q^-l,q^-2S-l,S+l]]q and [[q^ ,q^ -2S -2,6 + 2]]q. for 
< (5 < (7 — 1 by [9, Theorem 10]. The claim follows from 
Theorem 18. ■ 
The existence of the codes in i) are merely established by a 
non-constructive Gilbert- Varshamov type counting argument. 
However, the result is interesting, as it asserts that there exist 
for example [[6, 1, 1, 3]]g subsystem codes for all prime powers 
q > 7, [[7, 1, 2, 3]]q subsystem codes for all prime powers q > 
7, and other short subsystem codes that one should compare 
with a [[5,l,3]]q stabilizer code. If the syndrome calculation 
is simpler, then such subsystem codes could be of practical 
value. 

The subsystem codes given in ii)-vi) of the previous corol- 
lary are constructively established. The subsystem codes in 
ii) are derived from Reed-Muller codes, and in iii)-vi) from 



Subsystem Codes 



118,1, 5, 2]]3 
[[8,4,2,2]]3 
[[8,5,1,2]]3 

[[9,1,4,3]]3 

[[9,4,1,3]]3 



[[15.1, 10, 3JJ4 
[[15,9,2,3]]4 
[[15,10,1,3]]4 
[[16,1,9,4]]4 



[[24, 1,17,4JJ5 
[[24,16,2,4]]5 
[[24,17,1,4]]5 
[[24,19,1,3]]5 
[[24,21,1,2]]5 
[[23,1,18,3]]5 
[[23,16,3,3]]5 



[[48,1,37, 6jj7 



Paient 
Code (RS Code) 



[8,6,3]32 

[8.3.6) 32 

[8.2.7) 32 
[9,6,4]t 5 = 3 
[9,3,7U,,5 = 6 



[15,12,4J42 

[15.4. 12) 42 

[15.3.13) 42 
[16,12,5)^2, 5 = 4 



[24,20,5)52 

[24.5.20) 52 

[24.4.21) 52 

[24.3.22) 52 

[24.2.23) 52 
[23,20,4]*2,5 = 
[23,5, 19]%, S = 



-- 5 
20 



;42;7T 



72 



Punctured code 
t Extended code 



Reed-Solomon codes. There exists an overlap between the 
parameters given in ii) and in iv), but we list here both, since 
each code construction has its own merits. 

Remark 21. By Theorem 13, pure MDS subsystem codes can 
always be derived from MDS stabilizer codes, see Table III. 
Therefore, one can derive in fact all possible parameter sets 
of pure MDS subsystem codes with the help of Theorem 18. 

Remark 22. In the case of stabilizer codes, all MDS codes 
must be pure. For subsystem codes this is not true, as the 
[[9, 1,4, 3]]2 subsystem code shows. Finding such impure Fq- 
linear [[n, k, r, d]]q MDS subsystem codes with k + r = n — 
2d + 2 is a particularly interesting challenge. 

Recall that a pure subsystem code is called perfect if 
and only if it attains the Hamming bound with equality. 
We conclude this section with the following consequence of 
Theorem 18: 

Corollary 23. If there exists an Fq-linear pure [[n, k, d]]q 
stabilizer code that is perfect, then there exists a pure Fq- 
linear [[n, k — r,r, d]]q perfect subsystem code for all r in the 
range < r < k. 

VI. Extending and Shortening Subsystem Codes 

In Section IV, we showed how one can derive new subsys- 
tem codes from known ones by modifying the dimension of 
the subsystem and co-subsystem. In this section, we derive 
new subsystem codes from known ones by extending and 
shortening the length of the code. 

Theorem 24. If there exists an {{n, K, R, d))q Clifford sub- 
system code with K > 1, then there exists an ({n+1, K, R,> 
d))q subsystem code that is pure to 1. 

Proof: We first note that for any additive subcode X < 
Fq", we can define an additive code X' < Fq"+^ by 

X' ^ {{aa\bO) \ {a\b) eX,ae FJ. 
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We have \X'\ = q\X\. Furthermore, if {c\d) G X^^, then 
(cajdO) is contained in for all a in Fg, whence 

{X^^y C [X')^'. By comparing cardinalities we find that 
equality must hold; in other words, we have 

By Theorem 1, there are two additive codes C and D 
associated with an {{n, K, R, d))q Clifford subsystem code 
such that 

|C| = g"i?/A' 

and 

\D\ = \CnC^'^\=q"/{KR). 

We can derive from the code C two new additive codes of 
length 2n + 2 over F,, namely C and D' = C'n{C')^^. The 
codes C and D' determine a {{n + 1, K' , R' ,d'))q Clifford 
subsystem code. Since 

D' = c" n (C)-^- - c" n (C^O' 

we have \D'\ = q\D\. Furthermore, we have \C'\ = q\C\. It 
follows from Theorem 1 that 

(i) K' = q-+^/^/\CW\ = 'l''/^/\C\\D\ = K, 

(ii) R' = {\C'\/\D'\y/^ = i\c\/\D\y/^ = R, 

(iii) d' = swt((D')^^ \ C) > swi{{D^' \ Cy) = d. 
Since C contains a vector (Oa|00) of weight 1, the resulting 
subsystem code is pure to 1. ■ 

Corollary 25. If there exists an [[n, fc, r, subsystem code 
with k > and < r < k, then there exists an [[n+ 1, k, r, > 
d]]q subsystem code that is pure to 1. 

We can also shorten the length of a subsystem code in a 
simple way as shown in the following Theorem. 

Theorem 26. If a pure ((n, A', R, d))q subsystem code exists, 
then there exists a pure {{n — l,qK, R,d — subsystem 
code. 

Proof: By [^, Lemma 10], the existence of a pure Clifford 
subsystem code with parameters {{n, K, R,d))q implies the 
existence of a pure {{n,KR,d))q stabihzer code. It follows 
from [ 1 ' , Lemma 70] that there exist a pure ((n — 1, qKR, d — 
l))q Stabilizer code, which can be regarded as a pure {{n — 
1, qKR, l,d— l))q subsystem code. Thus, there exists a pure 
((?! — 1, qK, R, d—l))q subsystem code by Theorem 12, which 
proves the claim. ■ 
In bracket notation, the previous theorem states that the 
existence of a pure [[71, fc, r, subsystem code implies the 
existence of a pure [[n — l,k + l,r,d ~ l]]q subsystem code. 

VII. Combining Subsystem Codes 

In this section, we show how one can obtain a new subsys- 
tem code by combining two given subsystem codes in various 
ways. 



Theorem 27. If there exists a pure [[ni, ki,ri, c?i]]2 subsystem 
code and a pure [[ri2, ^21 ''21 c'2]]2 subsystem code such that 
^2 +''2 ^ '^i. then there exist subsystem codes with parameters 

[[ni +712- k2 - 7-2, /ci + ri - r,r, d]]2 

for all r in the range < r < ki + ri, where the minimum 
distance d > minjdi, di + d2 — fc2 ~ ^2}- 

Proof: Since there exist pure [[ni,ki,ri,di]]2 and 
[[n2, k2, r2, d2]]2 subsystem codes with k2+r2 < ni, it follows 
from Theorem 12 that there exist stabilizer codes with the 
parameters [[ni,fci + ri,di]]2 and [[n2,fc2 + r2,d2]]2 such 
that A;2 + r2 < rii. Therefore, there exists an [[rii + 712 — ^2 — 
r2, fci + Ti, d]]2 stabilizer code with minimum distance 

d > minjdi, di + c?2 ^ ^2 ^ ^2} 

by [ , , Theorem 8]. It follows from Theorem 10 that there 
exists [[ni +772 — ^2 — 7-2, /ci + ri — r,r,> d\]2 subsystem 
codes for all r in the range < r < fci + ri . ■ 

Theorem 28. Let Qi and Q2 be two pure subsystem codes 
with parameters [[n, fci, ri, dijjg and [[ti, fc2, '"2, ^2]]^, respec- 
tively. If Q2 ^ Qi, then there exists pure subsystem codes with 
parameters 

[[277, fci + + ri + r2 - r, r, d]]q 

for all r in the range < r < fci + fc2 + ri + 7*2, where the 
minimum distance d > min{di, 2^2}- 

Proof: By assumption, there exists a pure [[77, ki, r^, di]]q 
subsystem code, which implies the existence of a pure [[77, ki + 
ri,di]]q stabilizer code by Theorem 12, where i G {1,2}. 
By [11, Lemma 74], there exists a pure stabilizer code 
with parameters [[2r7, /ci + k2 + ri + 7'2,c?]]g such that d > 
mm{2d2,di}. By Theorem 10, there exist a pure subsystem 
code with parameters [[277, fci + fc2 + ri + 7-2 — r, r, d]]q for all 
r in the range < r < fci + fc2 + + '^2, which proves the 
claim. ■ 
Further analysis of propagation rules of subsystem code 
constructions, tables of upper and lower bounds, and short 
subsystem codes are presented in [2]. 

VIII. Conclusion and Discussion 

Subsystem codes are among the most versatile tools in 
quantum error-correction, since they allow one to combine 
the passive error-correction found in decoherence free sub- 
spaces and noiseless subsystems with the active error-control 
methods of quantum error-correcting codes. In this paper we 
demonstrate several methods of subsystem code constructions 
over binary and nonbinary fields. The subclass of Clifford 
subsystem codes that was studied in this paper is of particular 
interest because of the close connection to classical error- 
correcting codes. As Theorem 1 shows, one can derive from 
each additive code over F^ an Clifford subsystem code. This 
offers more flexibility than the slightly rigid framework of 
stabilizer codes. 

We showed that any F^-linear MDS stabilizer code yields 
a series of pure F^-hnear MDS subsystem codes. These 
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codes are known to be optimal among the -linear Clifford 
subsystem codes. We conjecture that the Singleton bound 
holds in general for subsystem codes. There is quite some 
evidence for this fact, as pure Clifford subsystem codes and 
Fq-linear Clifford subsystem codes are known to obey this 
bound. We have established a number of subsystem code 
constructions. In particular, we have shown how one can derive 
subsystem codes from stabilizer codes. In combination with 
the propagation rules that we have derived, one can easily 
create tables with the best known subsystem codes. Further 
propagation rules and examples of such tables are given in [2], 
and will appear in an expanded version of this paper 
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Appendix 

We recall that the Hermitian construction of stabilizer codes yields Fg-linear 
stabilizer codes, as can be seen from the following reformulation of [9, Corollary 2]. 

Lemma 29 ([')]). // there e.xists an F^2 -linear code X C F^2 ■s"'^'' ^ Q X''-'', 
then there e.xisl.^ an Fq-Unear code C C F^^ such that C C C'^^, \C\ — \X\, 
swt(C^= -C)= vit(X^h - X) and swt(C) = wt(X). 

Proof: Let {1,13} be a basis of F^2 /F,. Then tr^2/,(/3) = /3 -|- /S' is an 
element /3o of F,; hence, /3' = — /3 + /So- Let 

C = {{u\v) \u,v £F]^,u-\- pv £ X}. 

It follows from this definition that \X\ = \C\ and that wt(X) = swt(C). 
Fuithermore, if it -t- f3v and u' -t- f3v' are elements of X with u, v, u' , v' in F^ , then 

= (u + l3v)'' ■ (u' + 0v') 

— u ■ u' + ■ v' + /Sqv ■ u' + p(u ■ v' — V ■ u'). 

On the right hand side, all terms but the last are in F^; hence we must have {u • v' — 
V ■ u') — 0, which shows that {u\v) 1.^ (^'1^^'), whence C C C^^ . Expanding 
X^ti in the basis {1 /?} yields a code C' C , and we must have equality by a 

dimension argument. Since the basis expansion is isometric, it follows that 

swt(C-^= -C) = wt(X^'' - X). 

The Fg-linearity of C is a direct consequence of the definition of C . ■ 



